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Abstract 

Recent experiments have found that polyisocyanates composed of a mix- 
ture of opposite enantiomers fohow a chiral "majority rule:" the chiral order 
of the copolymer, measured by optical activity, is dominated by whichever 
enantiomer is in the majority. We explain this majority rule theoretically by 
mapping the random copolymer onto the random-field Ising model. Using this 
model, we predict the chiral order as a function of enantiomer concentration, 
in quantitative agreement with the experiments, and show how the sharpness 

of the majority-rule curve can be controlled. 

PACS numbers: 61.41. +e, 05.50.+q, 78.20.Ek, 82.90.+j 



1 



Typeset using REVT^ 



Cooperative chiral order plays a vital role in the self-assembly of ordered supramolecular 
structures in liquid crystals organic thin films and lipid membranes P-[T^. One 

particularly simple and well-controlled example of cooperative chiral order is in random 
copolymers. Recent experiments have found that polyisocyanates formed from a mixture of 
opposite enantiomers follow a chiral "majority rule" [0|. The chiral order of the copolymer, 



measured by optical activity, responds sharply to slight differences in the concentrations 
of the enantiomers, and is dominated by whichever enantiomer is in the majority. In this 
paper, we show that the majority rule can be understood through a mapping of the random 



copolymer onto the random-field Ising model ||T^-[T6|. Using this model, we predict the 
chiral order as a function of enantiomer concentration, in quantitative agreement with the 
experiments, and show that the sharpness of the majority-rule curve is determined by two 
energy scales associated with the chiral packing of monomers. 

In a series of experiments. Green et al. have investigated chiral order in polyiso- 
cyanates |T^. This polymer consists of a carbon- nitrogen backbone with a pendant group 
attached to each monomer, as shown in Fig. |I[ Although the backbone is nonchiral, steric 
constraints force the molecule to polymerize in a helical structure. If the pendant group is 
also nonchiral, the helix is randomly right- or left-handed. A long chain then consists of 
domains of fixed helicity, separated by occasional helix reversals. On average, there are equal 
right- and left-handed domains, leading to zero net optical activity. However, if the pendant 
group is chiral, there is a preference for one sense of the helix, which leads to a net optical 
activity. Because of the cooperative interaction between the monomers in a domain, even a 



very small chiral influence leads to a large optical activity [|T8[. Most recently. Green et al. 
have synthesized random copolymers with a mixture of right- and left-handed enantiomeric 
pendant groups, with concentrations p and 1 — p, respectively |13[ . The resulting optical 



activity, shown in Fig. ^, has a surprisingly sharp dependence on p. A 56/44 mixture of 
enantiomers has almost the same optical activity as a pure 100/0 homopolymer, and even a 
51/49 mixture has a third of that optical activity. 

To explain this cooperative chiral order theoretically, we map the random copolymer onto 



the one-dimensional random-field Ising model, a standard model in the theory of random 
magnetic systems |T^-|T6[. Although related models have been applied to other polymer 
systems |T9|-p2|, our theory gives a new, direct correspondence between the Ising order 



parameter and the optical activity. This correspondence provides a novel experimental test 
of predictions for the random-field Ising model. Let the Ising spin (Xj = ±1 represent the 
sense of the polymer helix at the monomer i. The energy of a polymer can then be written 



as 

N-l N 



H = - J aiai+i - hiai. (1) 



i=l i=l 



Here, the random field hi specifies the enantiomeric identity of the pendant group on 
monomer i: /ij = if it is right-handed (with probability p) and hi = — /i if it is left-handed 
(with probability 1 — p). This field is a quenched random variable; it is fixed by the poly- 
merization of each individual chain. The parameter 2h is the energy cost of a right-handed 
monomer in a left-handed helix, or vice versa. Molecular modeling gives 2h ^ 1.7 kJ/mol 
(0.4 kcal/mol) for the pendant group used in these experiments The parameter 2J is 
the energy cost of a helix reversal. Fits of the optical activity of pure homopolymers as a 
function of temperature T give 2J ^ 17 kJ/mol (4 kcal/mol) [|r^. Note that 2J is much 
greater than fc^T ^ 2.5 kJ/mol (0.6 kcal/mol), but 2h < ksT. The degree of polymerization 



ranges from 350 to 5800 in the experiments ||T3[. The magnetization of the Ising model, 

M^(lg.,). (2) 

corresponds to the chiral order parameter that is measured by optical activity. To predict 
the optical activity of the random copolymer as a function of enantiomer concentration, we 
must calculate M as a function of p. 

To calculate the order parameter M, we note that each chain consists of domains of 
uniform helicity cTj. As an approximation, suppose that each domain has length L, which 
is to be determined. Each domain responds to the total chiral field /itot = J2hi of the 
monomers in it. Because the domain is uniform, the response is M{htot) = tanh(/itot/^B7'), 



equivalent to a single spin in a magnetic field. Averaging over the probability distribution 
P{htot), we obtain 

M = 1" rf/itotP(/^tot) tanh (^^^ . (3) 
The probability distribution P{htot) is a binomial distribution. For large domains, it can be 



approximated by a Gaussian with mean 2hL{p — |) and standard deviation 2h[Lp{l — p)] 



1/2 



For p ~ |, the standard deviation becomes hL}^"^. Furthermore, if the width of the Gaussian 
is much greater than the width of the tanh, hL}^"^ ^ ksT , then the tanh can be approximated 
by a step function. The expression for M then becomes 



M ^ erf 



1/2 / 1 



(4) 



We must now estimate the domain size L. The domain size is determined by (a) the 
distance Lj-f between helix reversals that are induced by the random field, (b) the distance 
Lth between helix reversals that are induced by thermal fluctuations, and (c) the chain length 
A^. Because each of these effects contributes to the density 1/L of domain boundaries, we 
expect 

- ~ — — — (5) 

For p ~ ^ , the random-field domain size L^f can be estimated using a variation of the Imry- 
Ma argument for the random- field Ising model [jl4|. A domain forms when the field energy 



grows to equal the boundary energy J. By equating these two energies, we obtain the 
random-field domain size 



^rf 



(6) 



With the values 2 J ^ 17 kJ/mol (4 kcal/mol) and 2h ^ 1.7 kJ/mol (0.4 kcal/mol) appropri- 
ate for the experimental system, the random-field domain size becomes Lrf ~ 100 monomers. 
By comparison, the thermal domain size is Lth = e^'^^'^^^ ~ 800, and the chain length is 
A^ ~ 350-5800. Because Lrf is much less than Ltt and A^, we obtain L ^ Lj-f, i. e., the 
domain size is limited by random-field effects. 



To test this approximate calculation explicitly, and to obtain a more precise value of 
the domain size, we performed numerical simulations of the random-field Ising model. In 
these simulations, we used a series of chain lengths from 4 to 230, and used the values 
of J and h appropriate for the experiments. For each chain length, we constructed an 
explicit realization of the random field, then calculated the partition function and order 
parameter using transfer-matrix techniques. We then averaged the order parameter over 
at least 1000 realizations of the random field. Figure I shows M clS cL function of p for 
chain length = 230. These results can be fit very well to Eq. @, with the domain size 
L(230) = 96. The results for other values of N can be fit equally well. Figure ^ shows the 
fitted domain size L{N) as a function of chain length N . These results can be extrapolated 
using 1/L{N) = 1/Ljnax + 1/^, with L^ax = 164. This random- field domain size agrees well 
with Eq. (P), especially considering that the Imry-Ma argument is only a scaling argument. 
Thus, the chiral order parameter should indeed be given by Eq. (^, with the extrapolated 
domain size L = 164. 

To compare our theory with the experiment, we plot our prediction for the chiral order 
parameter M on top of the experimental data for the optical activity in Fig. |^. The prediction 
agrees very well with the data. In particular, M saturates at (p — |) ~ 0.06, a 56/44 
composition, in agreement with the data. This saturation point is a direct measure of 
1/(2L)^/^, which is controlled by the ratio of the two energy scales J and h. We emphasize 
that the theory matches the experimental data with no adjustable parameters, other than 
the relative scale of the optical-activity axis and the order-parameter axis. That relative 
scale is the optical activity of a pure 100/0 homopolymer. 

We can make two remarks about these results. First, both the random- field domain size 
Lrf and the thermal domain size Lth depend on p. The Imry-Ma argument above applies 
only to the regime where p ~ |. This is the appropriate regime for understanding the 
experiments, because all the significant variation in the chiral order parameter M occurs 
around p ~ Outside that regime, M saturates at ±1, and it is not sensitive to L. Second, 
the fact that Lrf <^ Lth shows that quenched disorder is much more significant than thermal 
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disorder in this copolymer system. Thus, the system is effectively in the low-temperature 
limit. If the temperature T were increased so that Lth ^ L^.^, then thermal disorder (i. e. 
entropy) would become more significant. However, that would require a temperature of 
T > 400 K, which is unrealistic because the system would degrade chemically. 

Using our theory, we can make two predictions for future experiments. First, the ex- 
periment could be repeated using different degrees of polymerization. For short chains, the 
domain size L is limited by the chain length A^, particularly for ^ 200, as shown in Fig. ^. 
Thus, shorter chains should give a broader version of the majority-rule curve. By contrast, 
longer chains should not give a sharper majority-rule curve, because the chains are already 
in the regime where L is approximately independent of A^. Second, the experiment could be 
repeated using pendant groups that are "less chiral," i. e. polyisocyanates with a lower en- 
ergy cost for a right-handed monomer in a left-handed helix. A lower value of the chiral field 
h should give a larger value of the domain size L ^ (J/h)"^, and hence a sharper majority- 
rule curve. (This prediction applies as long as Lrf <^ Lth, or 2h ^ 2Je^'^/'^^^ ~ 0.6 kJ/mol 
[0.14 kcal/mol]. Thus, h can be reduced by a factor of 3 from its value in the current ex- 
periments, and the majority-rule curve can become 3 times sharper. Beyond that point, the 
sharpness will be limited by the thermal domain size.) This second prediction might seem 
counter-intuitive, because one might expect a smaller chiral field to give a smaller effect. 
However, this prediction is reasonable, considering that the majority-rule curve is limited 
by the number of monomers that cooperate inside a single domain. If the local chiral field 
is reduced, then each monomer is more likely to have the same helicity as its neighbors, 
independent of the local chiral field, and hence the cooperativity increases. 

Finally, we note that the sharp majority-rule curve in polyisocyanates can be exploited 



in an optical switch If a mixture of enantiomers is exposed to a pulse of circularly po- 
larized light, one enantiomer is preferentially excited into a higher-energy state. That state 
can decay into either chiral form. Thus, a light pulse depletes the preferentially excited 
enantiomer and changes the enantiomer concentration p, which changes the optical activity. 
In other systems, this approach has been limited by the fact that light pulses induce only a 



slight change in p, and hence only a slight change in optical activity. However, in polyiso- 
cyanates, a slight change in p close to p = ^ is sufficient to induce a very significant change 
in optical activity. Indeed, this polymer has almost a binary response to changes in p, which 
is needed for an optical switch. Our theory shows how to optimize the majority-rule curve 
for use in an optical switch. 

In conclusion, we have shown that the cooperative chiral order in polyisocyanates can 
be understood through the random-field Ising model. The energy scales J and h, which 
arise from the chiral packing of the monomers, give the random-field domain size Lrf, which 
indicates how many monomers are correlated in a single domain. This domain size deter- 
mines the sharpness of the majority-rule curve. Our theory agrees well with the current 
experiments, and it shows how future experiments can control the chiral order. 

We thank M. M. Green and J. M. Schnur for many helpful discussions. 
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FIGURES 




FIG. 1. Molecular structure of the polyisocyanate with enantiomeric pendant groups derived 
from citronellic acid, which was studied in Ref. 



10 




FIG. 2. Symbols: Optical activity [a]^ of the random copolymer at the sodium D-line as a 



function of enantiomer concentration p, from Ref. |13|. Solid line: Theoretical prediction for the 
chiral order parameter M from Eq. with the domain size L = 164. The prediction agrees with 
the data with no adjustable parameters, other than the relative scale of the vertical axes. 
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FIG. 3. Numerical simulation of the chiral order parameter M as a function of enantiomer 
concentration p for chain length N = 230. The solid line shows a fit to the prediction of Eq. (Q), 
which gives the domain size L(230) = 96. 
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FIG. 4. The domain size L(N) from numerical simulations for several values of the chain 
length N. The solid line shows a fit to the extrapolation form 1/L{N) = l/Lmax + 1/-^, or 
L{N) = LjnaxN / {Liaax + N), which gives -Lmax = 164. The inset shows the extrapolation up to 
N = 6000. 
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